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The first-order bremsstrahlung emission spectrum is aduj/uj at zero temperature. If 
the radiation is emitted into a region that contains a thermal distribution of photons, 
then the rate is increased by a factor 1 -f N[uj) where N{u}) is the Bose-Einstein function. 
The stimulated emission changes the spectrum to aTdujjuP' for uj <^ T. If this were 
correct, an infinite amount of energy would be radiated in the low frequency modes. This 
unphysical result indicates a breakdown of perturbation theory. The paper computes the 
bremsstrahlung rate to all orders of perturbation theory, neglecting the recoil of the charged 
particle. When the perturbation series is summed, it has a different low-energy behavior. 
For iv ^ aT, the spectrum is independent of u and has a value proportional to du/aT . 
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I. INTRODUCTION AND SUMMARY 
A. Background 

In studies of the quark-gluon plasma to be produced in ultrarelativistic heavy ion coUi- 
sions and, more generally, in studies of field theory at finite temperature, a central concern 
is how the cancellation of infrared divergences affects various finite, physical quantities. 
The experimentally measured rates for certain low energy processes can be significantly 
modified by the infrared structure. 

This paper will investigate the emission of low energy photons by a charged particle 
that passes through a fixed-temperature plasma. The charged particle must undergo a 
collision in order to radiate. If the radiated energy is much smaller than the energy 
transfer in the collision then the inelastic cross section factors: 

^ da ^ dP(q^) da /,, , , 9,n 

^-dW"'"^5? (1.1) 

with the collision cross section da/dq^ independent of k ■ q, where q — p — p'. Thus to 
investigate the probability of radiation, P, the details of the hard scattering do not matter. 
All that matters is that the charged particle began with four-momentum and ended 
with four- momentum p'^. The radiation probability to first order in a is 



2uj 



dPi _ 
IFk ~ (27r) 

pol 



^j;|e^-jnMl + A^M] (1-2) 



where u = |^| and Niu) = [1 — exp(a;/T)] ^ is the Bose-Einstein function. When the 
radiated enegy is small, the matrix element of the electromagnetic current is 



where A is a momentum-space regulator that will be necessary later. The current re- 
sults from an on-shell charge of four-momentum p radiating an on-shell photon of four- 
momentum k, which gives a Feynman denominator (p — k)^ — m? = —2p-k. For low-energy 
radiation this current is valid regardless of the spin of the charged particle [1-3] . Because 
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of the denominators in (1.3) the radiation is mostly parallel to p or p'. When (1.2) is 
integrated over photon angles, the probability of radiating energy uj in any direction is 

^ = -[l + iV(a;)]e-/2^ (1.4) 

The charged particle momenta only occur in the function A: 

A{p-p') = -(-ln{\±^)-2] (1.5) 

TT \V 1 — V J 

where v is defined by p • p' = m^(l — v^)~^/^. In terms of momentum transfer Q, the 
limiting behaviors of A are [1,2] 

A^—^ for Q < m (1.6a) 

37r 

A ^ — In(^) for g » m (1.66) 

Formula (1.4) is classical, except for photon quantum statistics which produce the 
factor N . The formula fails at low energies because it predicts that the total energy 
radiated will be infinite: 

/■^-- dPi 

/ da;— — = oo (1- ' ) 

Jo cL; 

regardless of the value of -Emax- This does not occur at zero temperature, where the 
radiation probability is A/a;; but at non-zero temperature the Bose-Einstein factor makes 
it more singular: 

^ - f « T) (1.8) 

The breakdown comes from the smallest values of a;, where the classical current approxi- 
mation (1.3) works best. The resolution of the breakdown must occur from a higher order 
calculation. 
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B. Higher Orders 

To go to higher orders in a one must quantize the electromagnetic field with the 
interaction 

Li = -J^^{x)A^{x) (1.9) 

with the classical current (1.3) and the quantized field. This interaction will produce 
multiple emissions and absorptions of real photons as well as closed loops of virtual photons. 
For example, to second order in a there are three types of contributions to the probability 
of radiating energy lo: either radiate two photons whose energy totals or radiate one 
photon of energy uj + uj' and absorb another whose energy is a;', giving a net enegy of a;; or 
radiate a single photon of energy u with a one-loop correction. The infrared divergences 
cancel among these processes and give a finite answer dP2/du> to order a^. 

The semiclassical interaction (1.9) does not conserve energy-momentum: it allows 
the particle of momentum to radiate and still have momentum (provided the hard 
collision eventually defiects the charge to p'^) . This is a sensible approximation if each 
photon energy (real or virtual) is much smaller than the energy of the charged particle. 
To maintain this consistently in higher orders it is convenient to introduce the momentum 
cutoff A in (1.3) and require A -C (p^ -|- m^)^/^ and A -C (j)''^ + m^)^/^. 

Because the current is classical, one can compute the generating funtional for the 
multi-photon Green functions by performing a Gaussian functional integral. This allows 
one to calculate the exact multi-photon amplitudes /A. The probability of radiatiating a 
net energy uj to all orders in perturbation theory is 

— = J] / d^^...d^r.S{kU■■■k^-io) 

(1.10) 

lj]|>t(fci,...Ur 

pol 

Positive A;°'s correspond to emission of photons; negative A;°'s correspond to absorption of 
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photons. The phase space integration includes the appropriate statistical factors: 

2ki2nr \N{k) k^ = -k ^^-^^^ 

The amplitudes A4 contain infrared divergences from the virtual photon integrations, but 
the integration over the real photons in (1.10) produces an infrared-finite probabilty. This 
is the Bloch-Nordsieck cancellation at T ^ 0. The cancellation is very delicate: It would 
not occur if negative /c°'s (i.e. energy absorption from the heat bath) were omitted from 
the delta function in (1.10). 

C. Summary of Results 

Remarkably, one can calculate the probability (1.10) with only the approximation 

T -C A. The final answer obtained at the end of Sec. Ill is 

dP „M . a- e'"/='''e-l"l/'^ f2wT\^ 

dJ = '^<2+'2^)' 4.'TT(A) {—) (^-^2) 

Here ^ oc a is the same function oi p ■ p' as (1.5). This result also applies if u is negative, 
which means that the total emitted energy is less than the total absorbed energy. For 
u = —uj < 0, it has the property 

^ = e-^/^ ^ (1.13) 

In the first order result (1.4), this corresponds to exp(— a;/T)[l + N{Lij)] = N{uj). The zero 
temperature limit of (1.12) is discussed in Appendix A. 

At a; = the exact result (1.12) is finite. This resolves the total energy problem that 
arose in (1.7). The most interesting feature of (1.12) is that it involves two dimensionless 
quantities: A and u/T. To display a simpler form it is helpful to assume ^ <^ 1. This 
is quite reasonable because extremely large momentum transfers are necessary to make A 
large (for example, A = 1/16 requires Q/m = 10^ ). Even with A <^ 1 the behavior of 
(1.12) depends on the size of oj/T. If Att < uj/T then 

^ ~ ^[1 + Ar(a;)] e"'^/^ (AttT < a;) (1.14) 
duj ui 



which agrees with (1.4). However this does not hold at the smallest values of u. To 
approximate (1.12) for A -C 1 and u <^T one can use r(^) 1/z when \z\ <^ 1. Then 

HP AT 

Naturally (1.15) agrees with (1.14) in the region of overlap, AttT <^ lo <^ T. But for 
uj <^ AttT the radiation is suppressed relative to the first order rate. Fig. 1 compares 
(1.12) with the first-order result over a wide range of a;. 
Naturally the net radiated energy is finite, 

/ dw u;— ^ AT tan-\—) (1.16) 
Jo duj Att 

and numerically small. In the denominator of (1.15) the quantity AttT appears to be 
some type of thermal mass of order e^T. However, the folowing discussion will argue that 
AttT = r/2, where T is the radiation-damping rate. 

It is necessary to emphasize several points: (i) At T = the exact bremsstrahlung 
probability does not enjoy this damping. It is peculiar to T ^ 0. (ii) Whenever T > 
the suppression will occur for a charge of any energy. Whether the charge is relativistic or 
nonrelativistic alters the value of A, but the plateau shown in Fig. 1 will always exist. 

D. Interpretation of the Suppression 
One can understand the suppression of low frequency bremsstrahlung as the radia- 
tion reaction brought about by unitarity. First, we examine the failure of the first-order 
bremmstrahlung probability by writing it as 

du; 2tt , . 

1/2 (1-17) 



f,(.)=^(i±^) 



e l'^l/2A 



Let the Fourier transform of this be 

duj 



Then the bremsstrahlung probabihty has a time-dependence 

'"^^ = \m)\^ (1-19) 



dt 

As t — > oo the Fourier transform of (1.17) behaves hke 

Fi{t) -iVf (1.20) 

where F = 2tt AT. Therefore dPi/dt — > F as t — > oo. It is obviously unphysical for 
the radiation probabihty to remain constant after an infinitely long time. This disease is 
familiar in elementary atomic physics calculations of radiative transitions (e.g. 2p — > ls-l-7 
in hydrogen) . The resolution is that at higher orders the transition amplitude usually falls 
exponentially with time. Hence a reasonable guess is that higher order calculations should 
replace (1.20) by 

F{t) = -ie{t)Vf e-^*/2 (1.21) 
The Fourier transform to frequency gives 

UJ + iL /A 

The guess gives a radiation probability 

dP ? 1 



doj 27T + (F/2)2 



(1.23) 



which coincides with (1.15) the weak-coupling, low-energy limit. Note that this radiation 
damping does not come from modifications to the charged particle trajectory. 

The remainder of the paper describes the calculation that yields (1.12). Sec. II 
formulates the bremsstrahlung problem in detail and Sec. Ill performs the necessary 
integrations. Appendix A discusses the T = limit of (1.12) and its relation to conventional 
T = calculations. Appendix B discusses the relation of (1.12) to some previous work of 
mine. 



II. THERMAL BREMSSTRAHLUNG 

When charged particles are described by classical currents J'*, quantizing the radiation 
field becomes elementary. The generating functional for multi-photon Green functions is 

Z{J) = jD[A\exp[iJ d'^xiCo - J ■ A)] (2.1) 

where C is a contour in the complex plane that incorporates the temperature. For a 
real-time formulation [4-6] this gives 

ZiJ) = exp[^ J -0yJ^{k)D,,{k)r{-k)] (2.2) 

where the finite-temperature propagator is 

D^^k) = -g^, (^2^ - ^27r5(A;2)iV^ (2.3) 

with N = [exp(|/c|/T) - 1]-^ The amplitude for one real photon is 



M{k) = -ik\Jk)—--^ (27r)' 
= e^{k)J^{k) Z{J) 

This amplitude describes emission if k^ > 0; absorption, if k^ < 0. The amplitude for n 
photons is 

n 

M{ki,k2,...kn)^Z{J) l[[e^{ki)J>'{ke)] (2.5) 

To compute probabilities, one must square the amplitude and integrate over photon mo- 
menta. For the photon phase space it is convenient to employ the notation 

d^i^^6{kj) [e{k'i)+N{\k,\)] (2.6) 

This correctly weights photon emission with the statistical factor 1 + N; and photon 
absorption with the factor A^. The probability of radiating a net energy u> is 

n=l ^2.7) 

3 E 1-^(^1'-^^ 



pol 



The polarization sums give 



J2 \M{h, ...Ul' = \Z{J)f n Mke)J^{ke) (2.8) 

pol £=1 

It is important to note that the current (1.3) has the properties 

J^,{k)J^*{k) = J^ik)J^'i-k) 

(2.9) 

= -J^{k)J^{k) < 

Virtual Photons: The integration over all virtual photons is contained in the multi- 
plicative factor |Z(J)p. Using (2.2) and (2.3) gives 

|Z(J)|^ = exp(F) 

^ = ^1 2^'^'^(^)^'(^)[' + ''^] ('-'^^ 
Because J ~ 1/k, this integral contains both linear and logarithmic divergences in the 
infrared. 

Real Photons: The delta function that constrains the real photon energies can be 
represented as 

/OO J 
— e-^-^ (2.11) 
-OO 27r 

In (2.7) each integration over d'^kg will give the same function of z: 

R{z) = j d^e'^°^Y,M^)J^i^)\'' (2-12) 

pol 

More explicity this is 

d^k 



/d-^k 

X ([1 + Ar]e*'=^ + iVe-^'=^) 



(2.13) 



This is the contribution of the real photons. Because J ~ 1/k, this integral has both linear 
and logarithmic infrared divergence. 
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The probability (2.7) is 



■^e-"^^ exp[^(z)] 



— oo 



(2.14) 



(the n = term has no Fourier transform at u; 7^ 0) and R is 

R{z) = R{z) + V 

(2.15) 

X ([1 + N]e^^' + Ne-'^"" - [1 + 2N]) 

Each term in R has a clear interpretaion. Recall that z is the variable conjugate to the 
net energy uj. In (2.15) the term proportional to 1 + A?" represents the stimulated emission 
of energy; the term proportional to N represents the absorption of energy; the subtracted 
term 1 + 2 A'' represents all virtual photons (emitted and absorbed). At small k, J'^ 1/k 
and N 1/k. Expanding out the exp(±zA;2;) for small k shows that the linearly divergent 
terms dk/k^ cancel and the logarithmically divergent terms dk/k cancel. Thus i? is a 
completely finite function. (See (3.15) for R.) 



III. EXPLICIT INTEGRATION 
A. Angular Integration 
The remaining task is to perform the integrations necessary for the probability (2.14). 
The first step is to write 

^«=rf" (3.1) 

X ([1 + N]e'''' + e-'^' - [1 + 2Ar]) 
where A contains the angular integration: 

^= y/ -^J,{k)J^{k) (3.2) 
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Using k^^ = k{l, k) in the current (1.3) gives 

2{E'E-p'-p) 



{E'-p'.k){E-p.k) ^^^^^ 
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m m 



-k/A 



{E'-p'-kf {E-p-kf- 
The factors of k completely cancel so that A is independent of k. This is the same angular 
integration that arises in lowest order bremsstrahlung and is described in textbooks [2] . It 
turns out that A is a rather messy function of momentum transfer Q. It is a bit simpler 
to express it in terms of a = p' ■ p/m^: 



(3.4) 



As noted by Weinberg [7], it is simplest to express A in terms of the relative velocity 
V of the final particle in the rest frame of the initial particle (or vice versa) defined by 
p' ■ p = m'^{l - v^)~^/^ so that 

A{p.p') = -(-H\±^)-2) (3.5) 

B. Integration Over k 
To perform the k integration necessary for (3.1), separate R into a temperature- 
independent part and a temperature-dependent part: 

R{z) = {z) + Rt {z) (3.6a) 

Rq{z) = A / !^(e^fc^ - 1) exp(-A;/A) (3.66) 
Jo k 

^ / N , dk COS(kz) — 1 / , , . N /„ ^ X 

«-'^) = 2^i„ T exp(VT)-l '^"Pt-fc/A) (3.6c) 
Both Rq and Rt are infrared finite. For Rq{z), expand the integrand in powers of z 

Rq{z) =Ay ^ / dkk''-^ e-^^/^ (3.7) 
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The displayed integrals over k each give (n — 1)!A'^. The sum on n is elementary: 



R^{y) = -A\n[l - iAz] (3.8) 

For Rt{z), put k = 2Tx so that 

Rt{z) =A r ^ ^os(2r^a:)-l ^_,[i+, ^/^j 
Jo X smh(a;) 

This is finite at a; = but it is convenient to multiply the integrand by an additional 
convergence factor x^^ in order to use [8] 

/ da;x'^-i4T-- = 2i-'^r(^)Ck^] (3.10) 
Jo smh(a;) 2 

This gives 

Rt{z) = \l^ A2^-^V{^i) Re U[^iA] - (3.11) 

where q = q + iTz and q = 1 + T/ A. Now expand C ^ Taylor series about fx = using 
[8] 

C[;t/,Q]U=o = -g + l/2 (3.12a) 
^^l.=o = Mr(,)]-lln(27r) (3.126) 
When the limit ^ is taken in (3.11) the result is 

R{z) is the sum of (3.8) and (3.13). 

C. Integration Over z 
The last integral to perform is the Fourier transform 

dP r°° dz 



where the integrand is 



m( \\ - T{q + iTz)T{q-iTz) ^A 

ii-iKz)nq)nq) \ ^'-''^ 
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and q = 1 + T/A. Since A is non- integer, (3.15) has an infinite number of branch cuts in 
the variable z. For a; > (corresponding to net emission of energy by the charge) one can 
close the z contour in the lower half plane. Then there is a branch point at zq = —i/A 
and an infinte set of branch points at Zn = —i/A — in/T for n=l,2,3... that come from the 
poles of r(q' — iTz). Choose the branch cuts to run from z^ to 2;^ + oo as shown in Fig. 2 
so that 



n=0 ^" 

The discontinuity across the n'th branch cut is 



exp(i?(^)) 



DisCr 



exp(i?(2)) = [1 - e 



-i2^TA^ 



exp(i2(2)) 



Then put z = Zn-\- r where r is real: 



dP 



dw 



n=0 



dr 



^ ^e— ^exp(i?(2„ + r)) 



The value of the function along the branch cut is 



exp{R{zn + r)) ^ Cn{r) 
Cn{r) = 



^j, gi^(_n+l/2) 

Asinh(7rTr) 

r(l + 2T/A + n + iTr) 
r(l + n + zTr) [r(l + T/A)]2 

At this stage it is necessary to make the approximation T <^ A so that C„(r) 

allows the summation on n in (3.18) to be done: 

dP_ _ e'^/^^ e-'^/^sin(yl7r) (t^V" IM. 
U ~ sm{A7r/2 -iuj/2T) \~A^ J ^tT 



^0 



dre-''^''[smh{TTTr)\ 



This is a known integral [8] 



I{u) = 



r(l -^)r(^/2 + za;/27rT) 
7rT2i-^r(l - A/2 + iuj/2TTT) 
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(3.16) 



(3.17) 



(3.18) 



(3.19a) 

(3.196) 
1. This 

(3.20) 
(3.21) 



(3.22) 



Using the reflection property of the gamma function gives the final result quoted in (1.12): 

|r(o+^7^)l' (3-23) 



dw r(A)47r2T V ^ / 2 2tiT 
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APPENDIX: THE ZERO-TEMPERATURE LIMIT 
At zero temperature (1.12) has the simple behavior 

dP 1 /wx^ e"'^/^ 

(Al) 



duo uj \AJ T{A) 
and easily satisfies 



duj^ = l (A2) 
dw 



Since all detectors have some energy threshold for detection, physical probabilities must 
include an integration over the below-threshold photons. If the detector measures total 
radiant energy deposition (as in a calorimeter) above £^min then the appropriate probability 



IS 



n=l ■ 

This probability is 



Jo du; \ A J T{l + A) ^ ' 



with exp(— a;/A) neglected. 

The more familiar situation is a detector that is sensitive to single photons each of 
which has energy above some threshold £^min- In that case the single theta function in 
(A3) is replaced by a product: 

^(£^min - kl)d{E^in - /C2)...^(Smin - A^n) (^4) 
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This calculation is discussed by Weinberg [7] and gives a probability 

KA) {A6) 

where b{A) ^ 1 - (nA)^ /12 + ... is a different function that the 1/T{1 + A) that occurs m 
(A4). 

APPENDIX B: COMMENT ON A PREVIOUS CALCULATION 
In a previous paper on the semiclassical approximation [9], I computed the residual 
effects of infrared cancellation on processes which occur in a finite-volume heat bath but 
with no photons detected outside the heat bath. The absence of high energy photons 
means none were radiated by the charge. The absence of low energy photons, which have 
a short mean free path, is unavoidable because all those radiated by the charge would be 
thermalized in the heat bath and lose their identity. Consequently there is a threshold 
energy e that depends on the size of the system. In the notation of Sec. II, the photons 
with k < e contribute 

X ([1 + NY^^ + iVe"*^^) 
There is no constraint on the total energy u so the probability of no detected photons is 

/ duj^ = exp[i?,(0)] (S2) 

where 

RM = v + 



If e <^ T then -Re(O) —2AT/e. However, for reasonable mean free paths it invariably 
turns out that e > T, which made these thermal corrections small. 
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FIGURE CAPTIONS 

Fig. 1 The dimensionless probability, TdP/du, for a charged particle to radiate a net energy 
uj. The first order result fails when uj < AnT. For u> <^ AnT the exact result is 
constant, TdP/dw ^ I/Att'^. In this plot A = 0.01, which corresponds to a momentum 
transfer Q/m = 3. 

Fig. 2 Location of the branch cuts of the function exp(i?(2)) at z = —i/ A — in/T. For a; > 
the z integration contour may be wrapped around the branch branch cuts in the lower 
half-plane. 
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